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The Goldstone mode in the ordered phase of itinerant helimagnets, such as MnSi or FeGe, is 
determined and shown to have a strongly anisotropic dispersion relation. The softness of this mode 
is, in a well-defined sense, in between that of ferromagnetic and antiferromagnetic magnons, respec- 
tively. It is shown that this soft mode leads to nonanalytic corrections to Fermi-liquid behavior, 
with a T-contribution to the specific heat coefficient, and a T 5//2 -contribution to the resistivity. The 
quasi-particle inelastic lifetime shows anisotropic behavior in momentum space. 
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Helimagnets, which display long-ranged helical or spi- 
ral spin order, are less common, and less well known, then 
their more prevalent ferromagnetic and antiferromag- 
netic brethren. They nevertheless have received much 
attention lately, due to the unusual properties of the 
weak itinerant helimagnet MnSi. This material shows, 
at ambient pressure, helimagnetic order below a critical 
temperature T c ss 30 K . The wavelength of the helix 
is 2n/q w 180 A, with q the pitch wave number. Applica- 
tion of hydrostatic pressure p monotonically decreases T c 
until T c vanishes at p — p c w 14kbar 2\. In the ordered 
phase, neutron scattering shows the helical order with 
the helix pinned in the (11 Indirection due to crystal- 
field effects [3- In the disordered phase, quasi-static rem- 
nants of helical order are still observed at low tempera- 
tures close to the phase boundary, with an approximately 
isotropic distribution of the helical axis orientation |^|. 
In the entire disordered phase, up to a temperature of a 
few Kelvin, and up to the highest pressures investigated 
(~ 2p c ), pronounced non- Fermi- liquid behavior of the re- 
sistivity is observed, with the temperature dependence of 
the resistivity given by p(T — > 0) oc const. + T 3 / 2 |5(. In 
the ordered phase, on the other hand, the transport be- 
havior is Fermi-liquid-like, with a T 2 -dependence of the 
resistivity. No explanations have been given so far for 
these unusual properties. However, it is natural to spec- 
ulate that the remnants of helical order that are clearly 
observed in the paramagnetic phase have something to 
do with them. 

Given this situation, it is surprising that basic conse- 
quences of helical order, in particular the nature of the 
Goldstone mode that must result from the spontaneously 
broken symmetry, and its effects on observables, are not 
known. The present Letter addresses this issue. We will 
consider an isotropic itinerant helimagnet in its ordered 
phase. We will show that there is one Goldstone mode, 
the 'helimagnon', with a dispersion relation given by 



in the long- wavelength limit (|fe| < q). Here f2 is the fre- 
qency, k = (k±,k z ) is the wave vector, and q = (0,0, (?) 
is the pitch or axis vector of the helix, which we take 
to point in the z-direction. Notice that this dispersion 
relation is strongly anisotropic and much softer in the di- 
rection transverse to the pitch vector, Vt oc fe 2 for k z = 0, 
than in the z-direction, cx k z for k± = \k± \ = 0. 

Equation is our basic result, and all others follow 
from it. To the extent that the helimagnon is a well- 
defined quasi-particle, one expects its contribution to the 
internal energy to be U — £\ fi(fc) ra(f2(fe)/T). Here 
n{x) is the Bose distribution function, and we use unit 
such that h = k-Q = 1. It is easy to see that Eq. 
leads to a helimagnon contribution to the specific heat, 
C = dU/dT, given by C cx T 2 . Adding the leading 
Fermi-liquid contribution, which is linear in T, the low- 
temperature specific heat in the helical phase thus is 



C*(T->0) = 7 T + 7 2 T 2 + 0(T 3 In T). 



(2) 



fi(ft) cx Jk 2 z +k\/q 2 



(1) 



Here 7 is the usual specific- heat coefficient, and 72 de- 
pends on the properties of the helimagnon, see below. 

The single-particle Green function is non-diagonal in 
both spin and momentum space, but has two well-defined 
resonances that correspond to a split Fermi surface. The 
quasi-particle inelastic scattering rate, or inverse lifetime, 
due to scattering by helimagnons, l/r in (A;), one expects 
to be anisotropic due to the anisotropic nature of the 
helimagnon spectrum. Indeed, the poles of the Fermi 
surface at k = (0, 0, k z ) turn out to be 'hot spots' with a 
scattering rate l/r m cx T 2 , while elsewhere the scatter- 
ing is weaker with l/rj n cx T 5 / 2 . The transport lifetime, 
which determines the resistivity, averages over the Fermi 
surface, with the regions with a smaller scattering rate ef- 
fectively 'shunting out' or 'short circuiting' the hot spots 
0. As a result, the temperature dependence of the re- 
sistivity at low temperatures is given by 

p{T - 0) = Amp + 92 T 2 + P5/2 T 5/2 + 0(T 3 ). (3) 

Here p- lmv is the residual resistivity due to impurities, the 
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term proportional to T 2 is the usual Fermi-liquid contri- 
bution due to the screened Coulomb interaction between 
quasi-particles 0, and the T 5 / 2 contribution is due to 
scattering by helimagnons. Notice that the latter, albeit 
weaker than the leading Fermi-liquid term, is a nonana- 
lytic function of the temperature. 

We now sketch the derivation of these results; a com- 
plete account of the theory will be given elsewhere ||. 
Our starting point is a microscopic action So for free or 
band electrons, and an electron-electron interaction Si n t- 
The spin-singlet interaction will not be crucial for our 
purposes, and we suppress it for simplicity. The spin- 
triplet interaction has an isotropic, short-ranged con- 
tribution that can be modelled by the usual point-like 
Stoner interaction with an interaction amplitude T t 
In materials whose crystal structure lacks inversion sym- 
metry, the spin-orbit interaction leads to an additional 
chiral term [n], E3- This Dzyaloshinski-Moriya (DM) 
interaction is believed to be responsible for the helimag- 
netism observed in MnSi and FeGe ^3] • The spin-triplet 
interaction term can then be written 

f f 1/T 
s Lt = dxdy dr n l s (x, t) A io (x - y) n> s (y, r). 

(4a) 

Here n], (i = 1,2,3) are the components of the electronic 
spin density held which depend on position x and imag- 
inary time t, and 

Aij(x -y) = S i:j T t S(x-y) + e ijk C k (x - y), (4b) 

with eijk the Levi-Civita tensor and the C k the com- 
ponents of a vector-valued function of the position. In 
Fourier space, C k (p — > 0) = cp k - The coefficient c plays 
the role of a chiral coupling constant in the theory. 

A Hubbard-Stratonovich transformation that decou- 
ples the interaction term along the lines of Ref. Q leads 
to a Landau- Ginzburg- Wilson (LGW) functional for a 
helimagnet: 

5 hc n[M] = Sfm[M] +cj dxM(x) • [V x M(x)] , (5a) 

S im [M] = I J dxdy M(x)T(x~y)M(y) 

+ ljdx (M 2 (x)) 2 . (5b) 

The two-point vertex T reads, in Fourier space, 

rfotfin) = t + ap 2 + b\Q n \/\p\. (5c) 

Here we use a four-vector notation, x = (x, r) and J dx = 

f v dx Jq ■ Q n = denotes a bosonic Matsubara 

frequency, and M is the order parameter or Hubbard- 
Stratonovich field with components M{ (i = 1,2,3), 
whose expectation value is proportional to the magne- 
tization, t, a, 6, and u are the parameters of the LGW 



theory, they can be expressed in terms of the param- 
eters of the underlying microscopic theory if desirable. 
Sf m is Hertz's action for a quantum itinerant ferromag- 
net HEH. 

The action given by Eqs. 10 is easily seen to posses a 
helical saddle-point solution 



M sp (x) — m (cos(gz), shi{qz), 0) . 



(6) 



The free energy in saddle-point approximation is mini- 
mized by q = c/2a, and the amplitude m is determined 
by the saddle-point equation of state 



1 = cq ■ 



2T t 



V p tZ G o (P^ iLU n)G a {p-q,iu n 



Here A = mT t , and 



A 2 
(7) 



(8) 



is the inverse Green function for noninteracting electrons 



described by the action So- £ p 



ith €„ the 



energy-momentum relation and /i the chemical poten- 
tial, and uj n = 2nT(n + 1/2) is a fermionic Matsub- 
ara frequency. For explicit calculations we will use an 
effective-mass approximation, e p = p 2 /2m c , with m e the 
electronic effective mass. 

We now expand the action about the saddle-point so- 
lution, writing M(x) = M sp (x) + 8M(x). To Gaussian 
order in the fluctuations SM one obtains 

Sheii = S sp + ^ J dxdy SMi(x) % S(x - y) T t 

-e ijk S(x - y) r t c ^- - X : J cf (x, y) V 2 ] SM^y). (9) 

Here S sp is the saddle-point action, and 



(10) 



is a spin susceptibility for electrons in a reference ensem- 
ble described by an action for noninteracting electrons in 
an external field given by T t M sp (x). 

S rof = S - T t J dx M sp (x) ■ n s (x). (11) 

The Green function associated with this reference ensem- 
ble takes the form 

GkpiiUn) = S kp [<r+- a + (k,q;iuj n ) + a- + a-(k,q;iuj n )} 
+5k+ q , P v+b + (k,q;iLu n ) +8k-q, P °- b-(k,q;iu n ). (12a) 

Here a± = (<7i±i(72)/2 and a_| = er + (7_, a |_ = cr_cr + , 

with (7i 2 Pauli matrices, and 

a±(k,q;iv n ) = G ~ l (k±q,iLu n )/N±(k,q;iLu n ), (12b) 



b±(k,q;iw n ) = \/N±(k,q;iu n ), 



(12c) 
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FIG. 1: (a) The reference ensemble spin susceptibility, and 
(b) the lowest-order self-energy diagram contributing to the 
inelastic scattering rate. Solid and dashed line denotes the 
reference ensemble Green function and the helimagnon, re- 
spectively. 



with 

N±(k,q;iuj n ) = Go 1 (k,iaj n )Go 1 (k±q,iiu n )~X 2 . (12d) 

(. . .)r e f in Eq. (fTTTfl denotes a connected correlation func- 
tion with respect to SVef- Xrcf is given in terms of the 
reference ensemble Green function as shown in Fig. ^ 

The coefficients of the quadratic form given by Eq. 
© are thus known explicitly, and one can diagonalize it 
to find the excitations of the helically ordered state. We 
expect one soft or massless excitation, which corresponds 
to the Goldstone mode associated with the spontaneously 
broken symmetry in the ordered state. To find this mode, 
we parameterize the order parameter as follows |l4j 

(cos(q • x + 4>(x)) 
sin(q • x + 4>(x)) 
ipi(x) sin(q • x) + ip2{%) cos(q • x) 

' (13) 

and expand to first order in (j). The diagonalization prob- 
lem is now straightforward, although very cumbersome. 
The Goldstone mode, whose softness is ensured by the 
equation of state 10 is given by the linear combination 

g(k) =cj>(k,m n ) - i(k y /q)[l + 0(kl)}Mk) 

-i(W«)[l + O(fci)]lM*0, (14a) 

where k = (k,ifl n ). The helimagnon is the g-g correla- 
tion function, which for small frequencies (|f2 n | "C A) and 
small wave numbers (|fe| <C q) reads 



(g(k)g(-k)) = 



1 



1 



1 



2iV F A 2 2f(k) \f(k)+in n /2\ 
1 



/(fe) - in n /2\ 



where 



f(k) 



■kl/{2h 



(i 



k 4 J(2qk F ) 



(14b) 



(14c) 



with kp — y/2m c ^x the Fermi wave number and N-p = 
k F m e /2tt 2 . We see that the helimagnon is a propagating 
excitation with a dispersion relation given by Eq. (|T|). 

We now turn to the effect of this soft mode on various 
observables. The specific heat can be calculated either 
by the elementary method mentioned in the context of 



Eq. (PJ, or by calculating the free energy in saddle-point 
approximation. Either way one obtains 



9C(3) 3 1 T 2 
k F A 2 ' 



2V2n 



(15) 



where £ denotes the Riemann zeta function. 

The quasi-particle relaxation rate can be extracted 
from the imaginary part of the electronic self energy in 
the usual way. To linear order, the helimagnon contri- 
bution to the latter is given by the diagram shown in 
Fig. H For a quasiparticle on the Fermi surface, and for 
k = (fej_ = 0, k z ), one finds a scattering rate at asymp- 
totically low temperatures 



1 



1 _ 9V2fc F g 3 

r F (fe) ~ 64iV F m 2 

(2 + Q) 2 
+ [K z + (2 + Q) 2 ]3/ 2 



(2-Q) 2 



[K z + (2 



(16a) 



Here K z = \[2>k-pk z lm e \ and Q — \f%k-pq/2m e \. For 
k± the asymptotic scattering rate is smaller, 



1 9V2A kl-q 3 



1 



TF(fc) 



T 5/2 



4tt3 N F m 2 \k ± \VX Kfe+a+g) 4 

(16b) 

with A = 5.13 . . .. The crossover between these two types 
of behavior occurs at a temperature proportional to T x = 
(q/k F )(k 2 Jk 2 )X. 

To determine the temperature dependence of the resis- 
tivity we adapt a simple argument given for antiferromag- 
nets by Rosch in Ref. , who also showed that a varia- 
tional solution of the Boltzmann equation gives qualita- 
tively the same answer. Let To and To be microscopic 
temperature and time scales, respectively, t — T/Tq, 
and Ti mp the mean-free time due to impurity scattering. 
From Matthiessen's rule one has on the 'cold' parts of the 
Fermi surface, where the scattering rate goes like T 5 / 2 , 
l/fcoid = l/^imp + t 5 ^ 2 /T~Q. At the 'hot spots' one has 
1/Tiot = 1/Tmip + &t 2 /ro, with a a coefficient. From 
the dependence of the crossover temperature T x on k z 
and and applying elementary geometry (for small q 
the Fermi surface is approximately spherical) one finds 
that the effective area of the 'hot spots', normalized by 
the area of the Fermi surface, is linear in T, ahot — P t, 
with (3 another coefficient. The hot and cold areas of the 
Fermi surface are expected to act effectively like parallel 
resistors, i.e., the resistivity will be given by 



P = o i ahot Thot + (1 ~ ahot ) r coldJ 



n 



(17a) 



with n and e the electron density and charge, respec- 
tively. For the low-temperature behavior one thus finds 



P - Pimp (l + t 5/2 T imp /T + 0(t 3 ) 



(17b) 
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with p; mp = m c /ne 2 Ti mp the residual resistivity and p = 
m c /ne 2 To. This expression remains valid in the clean 
limit, when T- unp — > oo. Adding the usual Fermi- liquid 
contribution one obtains Eq. 1(3} . 

We finally discuss these results. The anisotropic nature 
of the dispersion relation, Eq. QJ, is analogous to that 
of the Goldstone mode in cholesteric liquid crystals [l6| . 
which also display helical order, although with a slightly 
different order parameter than a helimagnet. In cither 
case, the reason for the anisotropy is rotational invari- 
ance, which can be seen as follows 17]: A simple phase 
fluctuation, with tpi = ^2 = in Eq. (|13|l . would lead 
to a soft mode (the phase <j>) with fl 2 (k) oc k 2 . How- 
ever, an overall tilt of the helical axis, which cannot cost 
any energy, can be accomplished by choosing a phase 
4>(x) — 4>{x) — tpi x + cp2 V- To avoid this problem, the 
soft-mode action must not depend on Vj_</>, and the low- 
est order of k± allowed in the dispersion relation is k\. 

The helimagnon dispersion relation, Eqs. and (|14[l . 
is intermediate between isotropic ferromagnets (CI oc fe 2 ) 
and antiferromagnets (fi oc |fe|), respectively. In agree- 
ment with this observation, the contribution to the spe- 
cific heat is in between the one from ferromagnetic (C oc 
T 3 / 2 ) and antiferromagnetic (C oc T 3 ) magnons ^^|. Al- 
though it is weaker than the leading Fermi-liquid contri- 
bution, it is distinct from any contribution within Fermi- 
liquid theory, which does not produce a T 2 -term 0. 
Analogously, the helimagnon contribution to the tem- 
perature dependence of the resistivity is intermediate be- 
tween the one from ferromagnetic magnons (T 2 , Ref . H^l 
and antiferromagnetic magnons (T 5 , Ref. 20). Here the 
effect of the helimagnons is even more striking, as they 
produce a nonanalytic temperature dependence of the re- 
sistivity. Since the helimagnon contribution is stronger 
than the one from phonons, this effect should be observ- 
able if the T 2 background can be subtracted HU • 

All of these results are consistent with the experimental 
observation that the leading low-temperature behavior 
in the ordered phase of MnSi is consistent with Fermi- 
liquid theory We point out, however, that the result for 
the specific heat is more general and robust than the one 
for the resistivity. The former depends only on the exis- 
tence of well-defined quasi-particles with a dispersion re- 
lation given by Eq. (JIJ , whereas the latter corresponds to 
just a lowest-order variational solution of the Boltzmann 
equation. A more sophisticated transport theory, which 
takes into account either mode-mode coupling effects, or 
an interplay between quenched disorder and helimagnons 
beyond Matthiessen's rule, might well change the lead- 
ing temperature dependence of the resistivity. While the 
Boltzmann result seems to suggest that the non-Fermi- 
liquid behavior observed in the disordered phase is un- 
likely to be due to partial or quasi-static helical order, 
more work is needed to test this conclusion. 
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